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Abstract. We present generalisations of Wilson's theorem for 
double factorials, hyperfactorials, subfactorials and superfactorials. 



1. Introduction 

Wilson's theorem states that (p— 1)! = —1 Hp is prime, and {p—l)\ = 
otherwise, except for the one special case, p = 4. The result is 
attributed to John Wilson, a student of Waring, but it has apparently 
been known for over a thousand years; see [21], [SI Ch. II], [191 Ch. 11], 
[TOl Chap. 3] and [211 Chap. 3.5]. The first proof was given by Lagrange 
[16] . We give generalisations of Wilson's theorem for the four standard 
generalisations of the factorial function. 

Definition 1. For a natural number n, 

(a) the double factorial nil is the product of the natural numbers 
less than or equal to n that have the same parity as n, 

(b) the hyperfactorial H(n) is the number H(n) = YYk=i 

(c) the subfactorial \n is the number of permutations of the set 
{1, 2, . . . , r?,} that fix no element, 

(d) the superfactorial sf (n) is the number sf(?T,) = YYl^i k\ 

The idea underlying our approach is the following: if there are sensi- 
ble versions of Wilson's Theorem for these functions, then the numbers 
involved must be very special. If p is prime, then the obvious special 
numbers in the field Zp are 0, ±1 and when p is congruent to 1 mod 
4, the two square roots of —1. So the task is to verify that these are 
the values taken, and determine which value is taken in each case. We 
find that this simple approach for prime p works surprisingly well. For 
composite numbers, the hyperfactorial, subfactorial and superfactorial 
functions are all easily treated, while the double factorial holds some 
unexpected interesting surprises. We present the results for double fac- 
torials in Theorems |3l HI [6] and [71 The results for the superfactorial. 
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hyperfactorial and subfactorial are given in Theorems |2], [5] and M re- 
spectively. The following result is an unexpected consequence of this 
study. 

Theorem 1. If p is an odd prime, then modulo p 

{p - 1)!! = sf{p -1) = (-1)"^ H(p - 1). 

2. Factorials: double, hyper, sub and super 

According to the MacTutor History of Mathematics archive, the 
name factorial and the notation n\ were introduced by the French 
mathematician Christian Kramp in 1808 [12], but the symbol was not 
immediately universally adopted. In the English speaking world, the 
notation [n was still commonly used at the end of the 19*'^ Century 
[_5j . De Morgan wrote "Among the worst of barbarism^ is that of in- 
troducing symbols which are quite new in mathematical, but perfectly 
understood in common, language. Writers have borrowed from the 
Germans the abbreviation n! to signify 1.2.3 . . .{n — l).n, which gives 
their pages the appearance of expressing surprise and admiration that 
2,3,4,&c. should be found in mathematical results" |H]. 

Rev. W. Allen Whitwortlil apparently didn't share De Morgan's 
view. Whitworth introduced the subfactorial and a symbol for it, in a 
paper that begins with the words: "A new symbol in algebra is only 
half a benefit unless it has a new name. We believe that the symbol 
[n as an abbreviation of the continued product of the first n integers, 
was long in use before the name factorial n was adopted. But until it 
received its name it appealed only to the eye and not to the ear, and in 
reading aloud could only be described by a periphrasis" [25]. Subfac- 
torial n is the number of permutations of the set {1, 2, . . . , n} that fix 
no element. There are many symbols for the subfactorial. Whitworth 
used the symbol [n, in keeping with the notation for the factorial at the 



At some point, the word barbarisms was misspelt as barabarisms when quoted 
and this error has been reproduced in a great number of places. 

^ There is an amusing Australian connection concerning Whitworth. In the 
Cairns Post, 22 November 1890, an article reads: "The Rev. W. Allen Whitworth, 
in attempting a definition of gambling, as separable from mercantile speculation and 
legitimate enterprise, committed himself to the proposition that under ordinary 
circumstances and within the limits of moderation, it is (1) justifiable to back 
one's skill, (2) foolish to back one's luck, (3) immoral or fraudulent to back one's 
knowledge. There is no rule (4), but had there been it would doubtless have read 
- highly commendable to back Carbine". Carbine was the horse that won the 
Melbourne Cup in 1890. 
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time. These days n\ is sometimes used, but In seems more common. 
De Morgan would not be happy! 

The subfactorial has the following explicit formula. 

" (—IV 

(1) !^ = -!ES^- 

i=0 

The hyperfactorial function was also introduced in the 19^^ Century 
[TTj IH] and the superfactorial function followed at the turn of the 
Century [3], although the terms were introduced much later [23| [13]. 

Study of the double factorial goes back at least as far as 1948 [T7j, 
but it is probably much older. For even n, one has 

(2) n!!=2?-f! 

For n odd, n!! coincides with the Gauss factorial 712!, where in general 
n^! is the product of the natural numbers i < n that are relatively 
prime to m [7] . For a recent paper on the properties of double factorials 
and their applications, see P^ . 



3. Motivation 

Our motivation for considering the various factorial functions con- 
cerns the matrix 



A = 


n 

1 

1 


2 

22 

23 


3 .. 

32 .. 

33 .. 


■ P-l \ 
• {P - If 

■ ip-lf 




\i 


2P-1 


3f-^ . . 


■ {p-iy-'J 



There are many interesting facts and open questions related to the 
properties of the matrix A modulo p. For example, Giuga 's conjecture 
is that the sum of the entries in the bottom row is congruent to —1 if 
and only if p is prime [1] . 

Proposition 1. The above matrix A has determinant sf(p — 1). 
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Proof. Using the formula for the Vandermonde determinant [201 Chap. 1.2] 

we have 

/I 1 1 ... 1 \ 

12 3 ... p-1 
1 22 32 ... 



det(A) = {p-l)\ det 



(p 
{p 



\1 2P-2 3P-2 

n (^■-^) 

n o--i)!=sf(p-i). 



{p-iy-y 



2<j<p-l 



□ 



Remark 1. Notice that H(j9— 1) is precisely the product of the elements 
on the main diagonal of A. We will return to the matrix A briefly in 
Section [51 

4. The connection between the superfactorial and the 

double factorial 

Somewhat surprisingly, for prime p, the values sf(p — 1) and (p — 1)!! 
are congruent. Before proving this result, let us recall a well known fact 
which Lagrange proved in (TB] (see [121 H]). Let p be an odd prime. 
Then 

(3) (£^!)2^(_i)^ (inodp). 

Indeed, modulo p, one has p — i = —i for all i and so {p — 1)\ = 
(—1)^ (^i^O ^^"^ required claim follows from Wilson's Theorem. 

Theorem 2. If p is prime, then si{p — 1) = {p — 1)!! (mod p). 

Proof. The theorem is obvious for p = 2. In the following proof, p is 
odd and the congruences are taken modulo p. We have 

sf(p -l) = {p- 1)! {p - 2)! (p - 3)! . . . 3! 2! 1! 

(4) ={p- 1)!! ((p - 2)! (p - 4)! . . . 3! 1!)' . 
Note that for all 1 < z < p — 1 , 

_ -1 



(p- 

Hence 
(5) 



(p-l)(p-2)...(p-2 + l 
{p-i)\ = 



-l)(-2) 



1)! 
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So if is even, the factorials in @j cancel in pairs, giving sf(p — 1) = 
{p— 1)!!, as required. If is odd, cancellation of the factorials in @ 

leaves the middle term, giving sf — 1) = (p — 1)!! (^^!)^ and hence 
sf(p- 1) = (p- 1)!!, by dSD. □ 

Remark 2. In the above proof, the factorials in sf (p — 1) can also be 
cancelled in pairs using (E]), leaving only the middle one, and so 

(6) sf(p - 1) = (-l)SS'^ ■ = (-1)^ ■ 



Comparing this with (p— 1)!! = 2''2 -^y^-!, we obtain a simple derivation 



that 2'' 2 = (— 1)'^8 . In the following we will also require Euler's 

p— 1 

criterion Ch. Ill]: 2 is a quadratic residue if and only if = 1. 
Using this and Legendre's symbol we get a basic property of Gaussian 
reciprocity which allows to write: 

(7) (0 = (modp). 

5. The double factorial in the prime case 
If p is an odd prime, then (jS]) and give 

(8) ((p _ 1)11)2^ (£^1)2 (modp). 

In particular, {p — 1)!! = ±1 (mod p) when p = 3 (mod 4). The 
following result gives a more precise statement of this fact. 

Theorem 3. Suppose thatp is an odd prime andp = 3 (mod 4). Then 
{p—l)U = (— 1)*^ (mod p), where u denotes the number of nonquadratic 
residues j with 2 < j < | . 

Proof. Once again, the congruences will be taken modulo p unless oth- 
erwise stated. Let s := and Z+ := {1, 2, . . . ,p—l}. To evaluate the 
factor s! in ([7]) we use an argument that Mordell attributed to Dirichlet 
[T8] . Consider the involution (/p : x H- p—x on . Since p = 3 (mod 4), 
p cannot be expressed as a sum of two squares, and so 9? interchanges 
the set QR of quadratic residues with the set NR of nonquadratic 
residues. Let QR = {ri, r2, . . . , r^} and NR = {ni,n2, ■ ■ ■ , Ug}, with 
elements listed in natural order. Then 

s\ = rira . . . t^-at nirig . . . n^v = (-l)^rir2 . . . r^, 

where denotes the number of nonquadratic residues i < |. As p = 3 
(mod 4), one has rir2...rs = 1; see [22l p. 75]. Hence s! = (— 1)^. 
Thus by Remark [2], (p — 1)!! = 1 if and only if either 2 is quadratic 
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residue and is even, or 2 is nonquadratic residue and N is odd. Thus 
(p — 1)!! = 1 if and only if u is even. □ 

When p = 1 (mod 4), we have ((p — 1)!!)^ = —1 (mod p). In this 
case, because —1 is a quadratic residue, the involution ip used in the 
proof of Theorem [3] leaves the sets QR and NR invariant, so it is not 
useful. Instead we consider the involution -0 : x i— )■ of . It turns 
out that this approach is applicable for all odd primes. 

Theorem 4. Suppose thatp is an odd prime. Let /x denote the number 
of elements j less than | such that the inverse of j modulo p is also 
less than |. 

(a) If p = 3 (mod 4), then {p — = (modp). 

(b) If p = 1 (mod 4), then {p— 1)!! = (— (mod p), where ip 
is the unique natural number less than | with = —1 (mod p). 

Proof. We use the same notation as in the proof of Theorem [3J Notice 
that if J < I and j^^ < |, then the two terms cancel in ^y'-!. On 
the other hand, if j < | and j'^ > |, then p — < | and provided 
j P ~ ) the product i.{p — j^^) in gives —1. 

If p = 3 (mod 4), the number —1 is not a quadratic residue and so 
there is no j < | with j = p — In this case, b = = (—1)^, 

where s = Hence {p — 1)!! = ab = (— 1)^"^ ' ^ . Let p = 4A; + 3. 

Then 2^ + f = 2A;2 + 4A; + |. Hence {p - 1)!! = (-1)^ = (-1)"^, 
as required. 

If p = 1 (mod 4), we argue in the same manner, but now ip is the 
unique number with < | and ip = p — i'^ . Hence = (— 

where w = ^~2~^ . Thus {p — 1)!! = ab = ( — 1)^"8 ' i i^. Let 

p = 4:k + l. Then £^ + | = 2A;2 + 2k. Hence (p - 1)!! = {-l)^ip, as 
required. □ 

Together Theorems [3] and H] provide the following equivalence. 

Corollary 1. Ifp is an odd prime withp = 3 (mod 4), then the number 
V, of nonquadratic residues i with 2 < z < |, is even if and only if the 
number jj, of elements j less than | such that the inverse j^^ of j modulo 
p is also less than |, is congruent to 3 modulo 4- 

Remark 3. The results of Theorems [3] and |4] can be expressed in 
terms of class field numbers h, but the resulting statements are not 
as succinct as those given above. For the p = 3 (mod 4) case, one 
can use h{—p) = 2N — 1 (mod 4); see [18]. For p = 1 (mod 4), the 
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result can be expressed in terms of h{p) and the fundamental unit of 
the associated real quadratic number field; see [0]. 

We can now give the connection between the hyperfactorial and dou- 
ble factorial. 

Theorem 5. For p an odd prime, the hyperfactorial and double facto- 

p— 1 

rial are connected by the relation H(p— 1) = (— — 1)!! (mod p). 



Proof. Using Fermat's Little Theorem and Wilson's Theorem, we have 

{sf{p - 



p-i 

k 



H(p-i)=n^ 



fc=i 



sf - 2) 
[^i[p-\))v-\p-\) 

sf (p - 1) 
-1 -1 



sf(p-l) (p-1)!! 

by Theorem [21 By ([8]), we have 

(p— 1)!! : if p = 3 (mod 4) 
-(p— 1)!! : if p = 1 (mod 4). 



((p-1)!! 



,-1 



p— 1 

Thus H(p — 1) = (— l)~(p — 1)!!, as required. □ 

Remark 4. Note that by Proposition [H Remark[T]and Theorems [2] and 
[5l modulo p the determinant of the matrix A of Section [3] is (— 1)~ 
times the product of the elements on the main diagonal of A. 

6. Composite numbers 

When n is composite and n 7^ 4, one has (n — 1)! = (mod n). 
Similarly, it is obvious that sf(n — 1) =0 (mod n) and H(n — 1) = 
(mod n) for all composite natural numbers n. For the double factorial, 
the situation is more nuanced. The case of odd composites is not 
difficult. 

Theorem 6. If n is a composite odd natural number, then (n — 1)!! = 
(mod n) if n > 9, while 8!! = 6 (mod 9). 

Proof. Let n be a composite odd natural number. If n = ab, where 
a,b are co-prime, then a,b < so ^^! = (mod n). So we may 
assume that n is of the form n = p^ for some odd prime p, where 
k>2. If A; > 2, then Pjp''"^ are distinct and p,p''~^ < so ^! = 
(mod n). li n = p^ and p > 3, then p, 2p are distinct and p, 2p < 
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SO once again = (mod n). It remains to consider n = 9, which 
can be calculated by hand. □ 

However, when n is even, the pattern is less obvious. 

Theorem 7. Suppose that n = 2* (2 A; + 1), where i > 1 and k > 0. 
(a) ifi 



1, then (n — 1)! 

(b) ifi = 2, then (n- 1)! 

(c) ifi>2, then {n - 1)! 



= 2k + 1 (mod n), 
= -{2k + 1) (mod n), 



{2k 



(mod n) 



Proof. We first treat the case k = 0. Recall that for an arbitrary integer 
n, Gauss' generalisation of Wilson's Theorem [10^ Chap. Ill] (see also 
[2]) states that if / denotes the set of invertible elements in then 



n 

sei 



s = 



1 : if n = 4,p", 2p" {p an odd prime) 
1 : otherwise 



(mod n). 



For 77, = 2* the set of odd elements of Z„ is precisely the group of 
invertible elements in Z„. Hence Gauss' result gives 



[n 



1)! 



-1 



1 



if z = 2 
otherwise, 



(mod n) 



as required. 

Now assume k > 0. By the Chinese remainder theorem, the map 

: Zn ^ 1^2i X ^2fc+l 

m I-)- {ipi{m),ip2{m)) 

is a ring isomorphism, where ifi{m) (resp. <^2("^)) is the reduction of 
m modulo 2' (resp. 2k + 1). The inverse map is given by 

ip~'^{x, y) = 6?/2* + ax{2k + 1) (mod n). 

where a{2k + 1) + 62* = 1. The numbers a and b are defined modulo n. 
The proof focuses on the value of a. For i = 1 we may take a = 1. For 
i = 2, we can take a = {2k + 1). For i > 3, note that as the set of odd 
elements of is the group of invertible elements in 7j2^ ? so by Euler's 
Theorem, {2k + lf = 1 (mod 2*). Thus we may take a = (2/c + l)^' 
when i > 3. We have 



^{{n-m 



n 



X odd 



n 



\x odd 



((2*-l)!!,0) 
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Now from the k = case treated above 



(2*-l)!! = { ^ ■ ^ (modT). 

I 1 : otherwise. 

So when i = 1 we have 

(n - 1)!! = (p^\l, 0) = a{2k + 1) = {2k + 1). 

When i = 2 we have 

(n - 1)!! = ip-\-l,0) = -a{2k + 1) = -{2k + if. 

Finally, for i > 3 we have 

(n-1)!! = (/?-^(l,0) = a(2A; + l) = {2k+lf~\ 

7. The subfactorial 



□ 



For the (standard) factorial, the double factorial, the hyperfactorial 
and the superfactorial, the value at number n is congruent to zero mod 
n. It is in part for this reason that the value at n — 1 is of interest 
mod n. For the subfactorial however, the situation is different. Here 
the natural generalisation of Wilson's Theorem is the following. 

Theorem 8. If n is a natural number, then In = (—1)" (mod n). 

Proof. Modulo n we have from ([1]) 



J2'—^ = {-ir + n\J2—r^(-^r (modn) 



m = n\ 



□ 
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